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Abstract 

In terms of the stock exchange returns, we compute the analytic expres- 
sion of the probability distributions Fdax,+ and Fdax- of the normalized 
positive and negative DAX (Germany) index daily returns r(t). Further- 
more, we define the a re-scaled DAX daily index positive returns r(t) a and 
negative returns (— r(t)) a that we call, after normalization, the a positive 
fluctuations and a negative fluctuations. We use the Kolmogorov-Smirnov 
statistical test, as a method, to find the values of a that optimize the 
data collapse of the histogram of the a fluctuations with the Bramwell- 
Holdsworth-Pinton (BHP) probability density function. The optimal pa- 
rameters that we found are a + = 0.50 and aT = 0.48. Since the BHP 
probability density function appears in several other dissimilar phenom- 
ena, our results reveal universality in the stock exchange markets. 
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I. INTRODUCTION 



it is o 



The modeling of the time series of stock prices is a main issue in economics a nd finance anc 
a vit al imp ortance in t 



Lillo 



20011 1 and (IMantegna 



re ma nagement of large portfolios of stocks, see (IGabaix 



2001 



20011 ). Here we study de DAX indice. The DAX (Deutscher 



Aktien IndeX, formerly Deutscher Aktien-Index (German stock index)) is a blue chip stock 
market index that measures the development of the 30 largest and best-performing com- 
panies on the German equities market and represents around 80 % of the market capital 
authorized in Germany. The time series to investigate in our analysis is the DAX index 
from 1990 to 2009. Let Y(t) be the DAX index adjusted close value at day t. We define the 
DAX index daily return on day t by 



r{t) 



Y (t) - Y{t - 1) 
Y(t-l) 



We define the a re-scaled DAX daily index positive returns r(t) a , for r(t) > 0, that we call, 
after normalization, the a positive fluctuations. We define the a re-scaled DAX daily index 
negative returns (—r(t)) a , for r(t) < 0, that we call, after normalization, the a negative 
fluctuations. We analyze, separately, the a positive and a negative daily fluctuations that 
can have different s tatistical and economic natures du e, for instance, to t he leverage effects 



(see, for example, ([Andersen 



2004 



Barnhart 



20091 ) and flPintol . 120091 ) ). Our aim is to 



find the values of a that optimize the data collapse of the histogram of the a positive and 
a negative fluctuations to the universal, non-parametric, Bramwell-Holdsofworth-Pinton 
(BHP) probability density function. To do it, we apply the Kolmogorov-Smirnov statistic 
test to the null hypothesis claiming that the probability distribution of the a fluctuations is 
equal to the (BHP) distribution. We observe that the P values of the Kolmogorov-Smirnov 
test vary continuously with a. The highest P values P + = 0.19... and P~ = 0.24... of 
the Kolmogorov-Smirnov test are attained for the values a + = 0.50... and a~ = 0.48..., 
respectively, for the positive and negative fluctuations. Hence, the null hypothesis is not 
rejected for values of a in small neighborhoods of a + = 0.50... and a~ = 0.48.... Then, we 
show the data collapse of the histograms of the a + positive fluctuations and a~ negative 
fluctuations to the BHP pdf. Using this data collapse, we do a change of variable that 
allow us to compute the analytic expressions of the probability density functions fDAx,+ 



and Jdax,- of the normalized positive and negative DAX index daily returns 

//M.Y.. ■;•'•) = 5.58...X- ' 50 - W(22.2...x - 50 - - 1.99...). 
f D Ax,-(x) = 4.79...X" - 52 - W(20.12....x - 48 - - 2.01...) 

in terms of the BHP pdf Jbhp- We exhibit the data collapse of the histogram of the positive 
and negative returns to our proposed theoretical pdfs /dax,+ and Jdax-- Similar results are 
obser ved for some other stock indexes, prices of stocks, exchange rates and commodity prices 

cj)). Since the BHP pr obability density function appears in several 



see (IGoncalves 



2010b 



other dis similar phen omena (see, for instance, (jBramwelll . 



2002 



Dahlstedt 



2001 



2009b 



Pintd . 



Goncalvea . 



20101 )). our result reveals an universal feature of the stock exchange markets. 



II. POSITIVE DAX INDEX DAILY RETURNS 

Let T + be the set of all days t with positive returns, i.e. 

T+ = {t : r(t) > 0}. 

Let n + = 2524 be the cardinal of the set T + . The a re-scaled S&P100 daily index positive 
returns are the returns r(t) a with t G T + . Since the total number of observed days is 
n = 4758, we obtain that n + /n = 0.53. The mean /i+ = 0.09... of the a re-scaled DAX daily 
index positive returns is given by 

^ = 4 E r(tr (i) 

u ter+ 

The standard deviation cr+ = 0.045... of the a re-scaled DAX daily index positive returns is 
given by 



\ 



-1 Yl r(tr-" - «) 2 (2) 

71 t€T+ 



We define the a positive fluctuations by 

rUt) = ^4 (3) 

for every t G T + . Hence, the a positive fluctuations are the normalized a re-scaled DAX 
daily index positive returns. Let = —1.90... be the smallest a positive fluctuation, i.e. 

L t = mm{r+(t)}. 



Let = 5.51... be the largest a positive fluctuation, i.e. 

R+ = max{rj(t)}. 

We denote by F Qj+ the probability distribution of the a positive fluctuations. Let the truncated 
BHP probability distribution Fbhp,ci,+ be given by 

F BHPa+ (x) = F bhp(x) 

Fbhp(R£) — F B hp(L+) 

where Fbhp is the BHP probability distribution. 



We apply the Kolmogorov-Smirnov statistic test to the null hypothesis claiming that the 
probability distributions F a> + and F B HP,a,+ are equal. The Kolmogorov-Smirnov P value 
P a ,+ is plotted in Figured) Hence, we observe that a + = 0.50... is the point where the P 
value P a ,+ = 0.19... attains its maximum. 
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FIG. 1 The Kolmogorov-Smirnov P value P a ,+ for values of a in the range [0.4,0.6]. 

It is well-known that the Kolmogorov-Smirnov P value P a+ decreases with the distance 
||F a>+ - F B HP, a ,+ \\ between F a , + and F B HP,a,+- 

In Figure [2J we plot D a + + (x) = \F a + )+ (x) — F B HP,a+,+{ x )\ an d we observe that D a + + (x) 
attains its highest values for the a + positive fluctuations below or close to the mean of the 
probability distribution. 
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FIG. 2 The map A).50,+(» = 1^0.50,+ 0*0 ~ Fbhp,o.50,+ (x)\. 

In Figures [3] and HI we show the data collapse of the histogram f a +,+ of the a + positive 
fluctuations to the truncated BHP pdf fBHP,a+,+- 




FIG. 3 The histogram of the a + positive fluctuations with the truncated BHP pdf /bhp,o.50.+ on 
top, in the semi-log scale. 



Assume that t he probability dist ribution of the a + positive fluctuations r£+(t) is given by 
FBHP,a+,+, see ( Goncalves 
r(t) is given by 

a + x° 



fDAX,+ ( X 



2009al ). The pdf /dax,+ of the DAX daily index positive returns 



0"a+ ( F BHP ( R l+) - F BHP (fl 
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FIG. 4 The histogram of the a + positive fluctuations with the truncated BHP pdf /bhp,o.50.+ ° n 
top. 

Hence, taking a + = 0.50..., we get 

.//M.Y..i.r) = 5.58...X- ' 50 - W(22.2...x - 50 - - 1.99...). 

In Figures [5] and [6] , we show the data collapse of the histogram of the positive returns 
to our proposed theoretical pdf $dax,+- 
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FIG. 5 The histogram of the fluctuations of the positive returns with the pdf Jdax,+ on top, in 
the semi-log scale. 

III. NEGATIVE DAX INDEX DAILY RETURNS 

Let T~ be the set of all days t with negative returns, i.e. 

T~ = {t: r(t) < 0}. 
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0.02 0.04 



FIG. 6 The histogram of the fluctuations of the positive returns with the pdf /dax,+ on top. 

Let n~ = 2221 be the cardinal of the set T - . Since the total number of observed days is 
n = 4758, we obtain that n~ jn = 0.47. The a re-scaled DAX daily index negative returns 
are the returns (— r(t)) a with t G T~. We note that —r(t) is positive. The mean u,~ = 0.10... 
of the a re-scaled DAX daily index negative returns is given by 



— E (-^)) c 



n 



(4) 



t€T- 



The standard deviation a a = 0.050... of the a re-scaled DAX daily index negative returns 
is given by 

(5) 



\ 



— E (-^)) 2Q - (^) 2 



n 



teT- 



We define the a negative fluctuations by 

r~(t) = 



(6) 



for every t e T . Hence, the a negative fluctuations are the normalized a re-scaled DAX 
daily index negative returns. Let L~ = —1.94... be the smallest a negative fluctuation, i.e. 

K = min{r-(0}. 
Let R~ = 4.45... be the largest a negative fluctuation, i.e. 

R- = max{r-(t)}. 

We denote by F a _ the probability distribution of the a negative fluctuations. Let the trun- 
cated BHP probability distribution F BHP ^_ be given by 

Fbhp(x) 



BHP, 



a-(z) 



F~BHp(R a ) - F B Hp(L a ) 



where Fbhp is the BHP probability distribution. 

We apply the Kolmogorov-Smirnov statistic test to the null hypothesis claiming that the 
probability distributions F a _ and Fbhp,o> - are equal. 

The Kolmogorov-Smirnov P value P a ,- is plotted in Figure [7J Hence, we observe that 
a~ = 0.48... is the point where the P value P a - = 0.24... attains its maximum. 
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FIG. 7 The Kolmogorov-Smirnov P value P a - for values of a in the range [0.4,0.6]. 



The Kolmogorov-Smirnov P value P a _ decreases with the distance ||-F a ,_ — Fbhp,o,,-\\ be- 
tween F a _ and F B hp,<x,— 



In Figure [8], we plot D a -_(x) = \F a - ~(x) — F B HP,a--(%)\ and we observe that D a -_(x) 
attains its highest values for the a~ negative fluctuations above the mean of the probability 
distribution. 



In Figures 191 and [T0| we show the data collapse of the histogram / a -_ of the a negative 
fluctuations to the truncated BHP pdf fBHP, a -,-- 



Assume that th e probability distri bution of the a negative fluctuations r a _ (t) is given by 



Fbhp a- see (IGoncalvesl . 



2009al ). The pdf Jdax- of the DAX daily index (symmetric) 




FIG. 8 The map D 0A8 (x) = \F 0A8 (x) - F BHP0A8 (x)\. 



'"('0.48,-) 
■ ln l*RHPri4«J 




FIG. 9 The histogram of the a negative fluctuations with the truncated BHP pdf fBHP,0A8,~- on 
top, in the semi-log scale. 

negative returns —r(t), with T G T~, is given by 



fDAxA X ) 



a x 



°--(F BHP (R--)~F bhp (/£_))" 



Hence, taking a = 0.49..., we get 



f D AxA x ) = 1.7<)....r 11 ^ f mi r i2().\2 r"- 1 -- - 2.01...) 

In Figures [TT1 and [T2l we show the data collapse of the histogram of the negative returns 
to our proposed theoretical pdf foAx — 
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FIG. 10 The histogram of the a~ negative fluctuations with the truncated BHP pdf fBHP,o.A8,- 
on top. 



5 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 



FIG. 11 The histogram of the negative returns with the pdf /dax- on top, in the semi-log scale. 
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FIG. 12 The histogram of the negative returns with the pdf Jdax- on top. 



11 



IV. CONCLUSIONS 



We used the Kolmogorov-Smirnov statistical test to compare the histogram of the a pos- 
itive fluctuations and a negative fluctuations with the universal, non-parametric, Bramwell- 
Holdsworth-Pinton (BHP) probability distribution. We found that the parameters a + = 
0.50... and a~ = 0.48... for the positive and negative fluctuations, respectively, optimize the 
P value of the Kolmogorov-Smirnov test. We obtained that the respective P values of the 
Kolmogorov-Smirnov statistical test are P + = 0.19... and P~ = 0.23.... Hence, the null 
hypothesis was not rejected. The fact that a + is different from a~ can be do to leverage 
effects. We presented the data collapse of the corresponding fluctuations histograms to the 
BHP pdf. Furthermore, we computed the analytic expression of the probability distribu- 
tions Fdax,+ and Fdax- of the normalized DAX index daily positive and negative returns 
in terms of the BHP pdf. We showed the data collapse of the histogram of the positive and 
negative returns to our proposed theoretical pdfs fDAx,+ and jo ax-- The results obtained 
in daily returns also apply to other periodicities, such as weekly and monthly returns as well 
as intraday values. 



In (IPeixoto 



2001 



Goncalves 



2009a 



2010bl ). it is found the data collapses of the his- 



tograms of some other s tock indexes, price s of stocks, exchange rates and commodity prices 
to the BHP pdf and in (IGoncalvesl . l2010al ) for energy sources. 



Bramwell, Holdsworth and Pinton (IBramwelll . Il998l ) found the probability distribution 



of the fluctuations of the total magnetization, in the strong coupling (low temperature) 
regime, for a two-dimensional spin model (2dXY) using the spin wave approximation. From 



a statist ical physics point of view, one can thin 



system (jChowdhuryj . 



1999 



Gopikrishnanl . 



t that the stock prices form a non-equilibrium 



1998 



Lillo 



2001 



Plerou 



19991 ). Hence, the results 



presented here lead to a construction of a new qualitative and quantitative econophysics 
mod el for the stock ma rket based in the two-dimensional spin model (2dXY) at criticality 



(see (IGoncalvesl . 12010a )). 
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Appendix A: BRAMWELL-HOLDSWORTH-PINTON PROBABILITY DISTRI- 
BUTION 



The universal nonpa rametric BHP pdf was discovered by Bramwell, Holdsworth and 
Pinton (jBramwelll . 



19981 ). The BHP probability density function (pdf) is given by 



00 dx 



1 ^ 1 



1 v^iV-: 



JV-l 1 

3? 



.e 



EAT-l 
k=l 



77^7 -r — 4arctan ( 



(*)] ^-I'K 1 ^) 



(Al) 



where the {Afc}^ =1 are the eigenvalues, as determined in ( jBramwelll . l200ll ). of the adja- 
cency matrix. It follows, from the formula of the BHP pdf, that the asymptotic values 
for large deviations, below and above the mean, are exponential and double exponential, 
respectively (in this article, we use the approximation of the BHP pdf obtained by taking 
L = 10 and N = L 2 in equation ( 1Alj) ). As we can see, the BHP distribution does not 
have any parameter (except the mean that is normalize to and the standard deviation 
that is normalized to 1) and it is universal, in the sense that appears in several physical 
phenomena. For instance, the universal nonparametric BHP distribution is a good model to 
explain the fluctuations of order parameters in theoretical examples such as, mo dels of self- 
organ ized criticality, equilibrium critical behavio r, per colation phenom ena (see (jBramwelll . 
1998 )). the Snep pen model (see |Bramwell . 1998 ) and ( Dahlstedt . 2001 )). and auto- ignition 



fire models (see (pinha-Ray 



200ll )). The universal nonparametric BHP distribution is, also, 



an explanatory mo del for fluc t uatio ns of several phenomenon such as, width powe r in steady 



state systems (se e 



2001 



Goncalves 



2009b 



Bra mwelll . Il998l ) ) . fluctuations in river heights and flow (see (jBramwelll . 



c)), for the plasma density fluctuations and electrost atic turbulent 



fluxes measured at the scrape-off layer of the Alcator C-mod Tokamaks (see ( jVan Milligenl . 
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20051 )) and for Wolf's sunspot numbers fluctuations (see (jGoncalves 



2009dD ). 



REFERENCES 

Andersen, T.G., Bollerslev, T., Frederiksen, P. and Nielse, M., 2004, Continuos-Time Models, 
Realized Volatilities and Testable Distributional Implications for Daily Stock Returns Preprint. 

Barnhart, S. W. & Giannetti, A., 2009, Negative earnings, positive earnings and stock return 
predictability: An empirical examination of market timing Journal of Empirical Finance 16 
70-86. 

Bramwell, S.T., Holdsworth, P.C.W., & Pinton, J.F., 1998, Nature, 396 552-554. 

Bramwell, S.T., Fennell, T., Holdsworth, P.C.W.,& Portelli, B., 2002, Universal Fluctuations of the 

Danube Water Level: a Link with Turbulence, Criticality and Company Growth, Europhysics 

Letters 57 310. 

Bramwell, S.T., Fortin, J.Y., Holdsworth, P.C.W., Peysson, S., Pinton, J.F., Portelli, B. & Sellitto, 
M., 2001, Magnetic Fluctuations in the classical XY model: the origin of an exponential tail in 
a complex system, Phys. Rev E 63 041106. 

Chowdhury, D. and Stauffer, D., 1991, A generalized spin model of financial markets Eur. Phys. 
J. B8 477-482. 

Dahlstedt, K., & Jensen, H.J., 2001, Universal fluctuations and extreme-value statistics, J. Phys. 
A: Math. Gen. 34 11193-11200. 

Dynamics, Games and Science., 2010, Eds: M. Peixoto, A. A. Pinto and D. A. Rand. Proceedings 
in Mathematics series, Springer- Ver lag. 

Gabaix, X., Parameswaran, G., Plerou, V. & Stanley, E., 2003, A theory of power-law distributions 
in financial markets Nature 423 267-270. 

Gongalves, R., Ferreira, H. and Pinto, A. A., 2009a, Universality in the Stock Exchange Market 
Journal of Difference Equations and Applications (accepted). 

Gongalves, R., Ferreira, H. and Pinto, A. A., 2010a, Universality in energy sources, IAEE (Inter- 
national Association for Energy economics) International Conference (accepted). 

Gongalves, R., Ferreira, H. and Pinto, A. A., 2010b, Universal fluctuations of the Dow Jones 
(submitted). 

14 



Gongalves, R., Ferreira, H. and Pinto, A. A., 2010c, A qualitative and quantitative Econophysics 

stock market model (submitted). 
Gongalves, R., Ferreira, H., Pinto, A. A. and Stollenwerk, N., 2009b, Universality in nonlinear 

prediction of complex systems. Special issue in honor of Saber Elaydi. Journal of Difference 

Equations and Applications 15, Issue 11 &: 12, 1067-1076. 
Gongalves, R., and Pinto, A. A., 2009c, Negro and Danube are mirror rivers. Special issue Dynamics 

& Applications in honor of Mauricio Peixoto and David Rand. Journal of Difference Equations 

and Applications. 

Gongalves, R. Pinto, A. A., Stollenwerk, N., 2009d, Cycles and universality in sunspot numbers 

fluctuations The Astrophysical Journal 691 1583-1586. 
Gopikrishnan, P., Meyer, M., Amaral, L. & Stanley, H., 1998, Inverse cubic law for the distribution 

of stock price variation The European Physical Journal B 3 139-140. 
Lillo, F. and Mantegna, R., 2001, Ensemble Properties of securities traded in the Nasdaq market 

Physica A 299 161-167 (2001). 
Mantegna, R. & Stanley, E., 2001, Scaling behaviour in the dynamics of a economic index Nature 

376 46-49. 

Pinto, A. A., 2010, Game theory and Duopoly Models Interdisciplinary Applied Mathematics, 
Springer- Ver lag. 

Pinto, A. A., Rand, D. A. and Ferreira, F., 2009, Fine Structures of Hyperbolic Diffeomorphisms. 

Springer- Verlag Monograph. 
Plerou, V., Amaral, L., Gopikrishnan, P., Meyer, M. k, Stanley, E., 1999, Universal and Nonuni- 

versal Properties of Cross Correlations in Financial Time Series Physical Review Letters 83 7 

1471-1474. 

Sinha-Ray, P., Borda de Agua, L. & Jensen, H.J., 2001, Threshold dynamics, multifractality and 
universal fluctuations in the SOC forest fire: facets of an auto-ignition model Physica D 157, 
186-196. 

Van Milligen, B. Ph., Sanchez, R., Carreras, B. A., Lynch, V. E., LaBombard, B., Pedrosa, M. A., 
Hidalgo, C, Gongalves, B. & Balbm, R., 2005, Physics of plasmas 12 05207. 



15 



